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ABSTRACT 



Synchrotron radiation mechanism, when electrons are accelerated in a relativistic shock, is 
known to have serious problems to explain the observed gamma-ray spectrum below the peak 
for most Gamma-Ray Bursts (GRBs); the synchrotron spectrum below the peak is much softer 
than observed spectra. Recently, the possibility that electrons responsible for the radiation cool 
via Inverse Compton, but in the Klein-Nishina regime, has been proposed as a solution to this 
[*T"| ■ problem. We provide an analytical study of this effect and show that it leads to a hardening 

of the low energy spectrum but not by enough to make it consistent with the observed spectra 
for most GRBs (this is assuming that electrons are injected continuously over a time scale 
comparable to the dynamical time scale, as is expected for internal shocks of GRBs). In par- 
ticular, we find that it is not possible to obtain a spectrum with a > —0.1 (/„ oc v a ) whereas 
the typical observed value is a ~ 0. Moreover, extreme values for a number of parameters 
are required in order that a ~ —0.1: the energy fraction in magnetic field needs to be less 
than about 10~ 4 , the thermal Lorentz factor of electrons should be larger than 10 6 , and the 
radius where gamma-rays are produced should be not too far away from the deceleration ra- 
dius. These difficulties suggest that the synchrotron radiation mechanism in internal shocks 
does not provide a self-consistent solution when a > — 0.2. 
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1 INTRODUCTION 

The dissipation mechanism responsible for the prompt emission of Gamma-Ray Bursts (GRBs) remains unknown. There have been various 
ideas put forth to explain it (for a review, see Piran 1999, 2004; Meszaros 2006; Gehrels et al. 2009). One of the main problems is to explain 
^ ■ the fact that the majority of GRBs exhibit a spectrum /„ oc v a , with a ~ below the peak of the spectrum (Preece et al. 2000), whereas 
the simplest version of the synchrotron model (in the so-called "fast cooling regime") predicts a — —1/2 (see, e.g., Ghisellini et al. 2000). 
Recently, a modified version of the synchrotron model, in which electrons that radiate below the peak of the spectrum cool via Inverse 
Compton (IC) in the Klein-Nishina (KN) regime (Derishev et al. 2003) has gained popularity (Bosnjak et al. 2009, Nakar et al. 2009, Wang 
et al. 2009, Fan 2010, Daigne et al. 201 1). The main idea is very simple: Electrons cooling via synchrotron mechanism (or IC in the Thomson 
regime) exhibit an energy loss rate oc , where 7 e is the electron Lorentz Factor (LF) and 8 = 2. The observed synchrotron spectrum is 
then f v . However, when the cooling of electrons is dominated by the IC in the KN regime, where the photon-electron 

interaction cross section scales as ~ 7" 1 , then 5 ~ 1 and /„ ~ v°. In this paper, we investigate this scenario in detail and explore its 
consequences. 

We analytically study the IC cooling in the KN regime assuming that electrons are injected continuously over a time scale comparable to 
the dynamical time scale, as is expected for the internal shock model of GRBs (Piran, Shemi & Narayan 1993; Katz 1994; Rees & Meszaros 
1994). Recently, Daigne et al. (2011) have provided a detailed numerical calculation of the same context; however, they have assumed that 
electrons are injected instantaneously in the internal shock. Therefore, our work and the work of Daigne et al. (2011) are complementary. 
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Essentially, Daigne et al. (201 1) deals with the case when electrons are no longer being injected and they simply cool, which happens when 
the internal shock has already passed through the shell. They consider the superposition of emission of many shells, for which the shock has 
crossed all of them, and the shells simply adiabatically expand and cool. We, however, consider the shock as it traverses the shell, accelerates 
electrons, and these radiate. 

The scenario presented in this paper has been considered before (Nakar et al. 2009; Fan 2010) and our results are consistent. However, in 
contrast with these works and with the work of Daigne et al. (201 1), the work presented here is applied to the prompt phase data of particular 
GRBs with a ~ 0: we analyze the data of GRB 080916C, a burst detected by the Fermi Satellite, in the context of the scenario described 
above, and provide constraints on this scenario based on available >100 MeV data. 

Recent developments on prompt GRB theory have cast doubt on the internal shock model (see, e.g., Kumar & Narayan 2009, Zou et 
al. 2009). New alternative models have been proposed to solve the low-energy spectral index problem described above and other prompt 
theory issues (see, e.g., Meszaros & Rees 2000, Drenkhahn & Spruit 2002; Lyutikov & Blandford 2003; Giannios 2008; Narayan & Kumar 
2009; Kumar & Narayan 2009; Lazar et al. 2009; Beloborodov 2010; Vurm et al. 201 1 ; Meszaros & Rees 201 1; Ioka 2010; Ioka et al. 201 1; 
Zhang & Yan 201 1, Bosnjak & Kumar 2012; Pe'er et al. 2012). It is, however, still relevant to critically test the internal shock model, in the 
particular case where electrons cool via IC in the KN regime, to assess its feasibility. 

It is important to mention that there exists a fraction of GRBs with < a < 1/3 (only about 25 per cent of GRBs have more than 50 
per cent of their spectra with a in this range; see Kaneko et al. 2006), that is, with spectra consistent with synchrotron radiation mechanism; 
however the scenario presented here is unable to explain them. In this work we focus on the majority of GRBs, which have a ~ 0; in 
particular, we study the case of GRB 080916C, which shows a = -0.02 ± 0.02 for most of its duration (Abdo et al. 2009). 

We set up our model and present the relevant time scales in Section 2. In Sections 3 and 4, we calculate the effect of IC cooling in the 
KN regime on the electron energy distribution and on the observed spectral slope, respectively. In Section 5 we derive the relevant physical 
parameters (radius of emission and total luminosity). In Section 6 we apply our results to GRB080916C, and to an average long-duration 
GRB. In Sections 7 and 8, we present a Discussion and our Conclusions. 



2 ELECTRON COOLING 



in observer frame is u v . We define 



Let us consider a GRB jet that has bulk LF V and bolometric 7-ray luminosity (isotropic equivalent) L. The peak of the GRB spectrum (uf v ) 

(1) 

where z is the redshift and h, m e and c are Planck's constant, the electron mass and the speed of light, respectively. 

The case we are trying to explain using the idea of electron cooling via IC cooling in the KN regime is when the observed spectrum 
below the peak is /„ oc v a , with a ~ 0; a is known as the low energy spectral index. We take a to extend from at least 10 keV to v v \ 10 
keV is roughly the lower energy limit of the GBM detector on board the Fermi satellite. 

Let us take the thermal LF of electrons that produce 10 keV photons (via synchrotron process) to be 74, and the LF of electrons producing 
photons of frequency v v to be 7i. We start with the assumption that electrons with LF 74 < 7 e < 7; cool primarily via IC in the KN regime. 
In order to satisfy this assumption, 74 should be such that 

74(l + ^ P >mec2) (2) 

which leads to the condition that 74 > £p Y. We define a variable 774 (7/4 ^ 1) which tells us how deep electrons of 74 are in the KN regime. 
With this, 74 is 



74 = T]4£p 



'r. (3) 



We can find the magnetic field strength in the jet comoving frame, B, so that electrons with LF 74 have synchrotron radiation at 10 keV 
in the observer frame. The observed synchrotron frequency of electrons of 74 is 

eB 7 |r „ ir » tr , B 7 |r 



= (1.15 x 10 eV) '* , , (4) 
27rm e c(l + 2) ' (1 + z) 

where e is the electron charge. For v = 10 keV, and using {3), we find the magnetic field 

B = (8.7 x 10 2 G)(1 +z)e^r 2 rj 3 , (5) 

where here and throughout the paper we use the usual notation Q n = Q/10™, with the exception of 7 e , r\ and the Compton-y parameter, Y; 
in these cases the subscript indicates the logio of the observed synchrotron frequency in eV we are referring to. 
The cooling time due to synchrotron radiation for an electron of LF y e , in the jet comoving frame, is given by 
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where ot is the Thomson cross section. For 7 e = 74 we find, using (O and (O, that 

t'syn = ( Lls )( 1 + Z j , 2 e 3- ^ 

We now calculate the electron cooling time due to IC scattering of 7-ray photons. The cross section for scattering photons of frequency 
Up by electrons of 74 is smaller than the Thomson cross section by a factor of w 774. Thus, 

74((TtL/774) 74^53 

or by substituting ((3), we find 

t' IC a (0.15s)R? 5 r 3e pL 5 - 3 1 . (9) 

Note that the IC cooling time is essentially independent of j e in the KN regime. The reason is that the energy of electrons is rn e c 2 y e and 
the IC power in the KN regime is approximately oc 7 e , therefore, the time scale is almost independent of y e . We will compare these cooling 
time scales with the dynamical time in the jet comoving frame 

t' dyn = = (33s)R 15 r 3 - 1 . (10) 



3 EFFECT OF IC COOLING IN KN REGIME ON ELECTRON DISTRIBUTION 

The electron energy distribution, n e , in steady state, for electrons of LF "/ e , is determined from the continuity equation 

A [7en( 7e )] = 5( 7 e) oc ( ffl " 7e J 7 < , (11) 

Ole I 7 e < 7i, 

and the cooling of electrons is determined by 

_ . _ arB 2 ^ + oknLjI , 12 - 
e 67rm e c 47ri? 2 r 2 m e c 2 ' 



where okn is the KN cross section which we write as <tkn = ct/ (77), and 



1 + 77 f 27,(1 + 77) ln(1 + 2?7 ^ , ln(l + 277) 1 + 37, 



(13) 



I + 277 / 2?, (1 + 2t?) 2 

We also define 7, as 77 = e P 7 e /r analogous to ®, and it indicates how deep electrons of 7 e are in the KN regime, defined in In this 
section we keep the dependence of 77 on j e . As a reminder, 77 was defined for a specific 7 e = 74 in the previous section; we will return to 
that same definition later on. With the use of d), equation can be rewritten as 

-%=7^- + ^-m, (14) 

where Tiyu an d T' IC are defined as 

T' = 6%meC (15) 
yn ~ a T B 2 ' 

and 

, _ 4nR 2 T 2 m e c 2 

J-IC = r ■ ( - 16 - ) 

OT Li 

Note that t' syn = T' ayn / r y e is the synchrotron cooling time, defined in ((6}, and t' IC KN = T' ic ot jipKN^e) = 2/c/(/( r ?)7e) is the IC 
cooling time in the KN regime (see eq. {8} and note that 7(77) ~ 77" 1 for 77 ^> 1). 

We can define Ykn = ^ ay nli'ic,KN = T' syn / {T' IC / /(»?)), and identify it as the Compton-F parameter in the KN regime for electrons 
with LF 7 e , which are deep in the KN regime as characterized by their parameter 77. The electron cooling rate is now 

-% = t^-(1 + Ykn). (17) 

J syn 

Since we are interested in the flux below the peak, we consider electrons with LF 7 e < 7*. Moreover, we need to consider the case 
7c < 7e < 7i, where 7 C is the cooling LF, which is the LF of electrons that cool on a dynamical time, td yn - The case 7 e < 7i < 7c gives 
fu oc v 1 / 3 below the peak, which is impossible to obtain for electrons accelerated in shocks for GRBs (Ghisellini et al. 2000; Kumar & 
McMahon 2008). 

The solution of the continuity equation dl lb for the the electron distribution for 7 e < 7, is n e ccj e ~ l and, using dl7t . the corresponding 
power-law index of the distribution is 
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Pi 



dlnn e 



dln7 e 



= 2 + 



Y KN d ^l 



f(v)(i + Y KN y 



As mentioned before, 77 depends on 7 e , and using equation d 1 3 fc we find 



drj 



2r]{2rf- - 39T/ 3 - 63t? 2 - 34?? - 6) - (2rj + l) 3 (rj 2 —477 — 6) 1h(2tj + 1) 



(18) 



(19) 



2»7 4 (2?7 + l) 3 

For 77 1, which we consider here, eq. ( I19t always yields a negative value, and since Ykn ^ 0, then 1 ^ pi ^ 2. Since we are 
interested in f„ oc u°, we will focus on finding the corresponding 77 and Ykn that give pi — ¥ 1, so that a — ¥ 0. Let us first calculate a for a 
given pi . 



4 SPECTRAL SLOPE WHEN ELECTRON DISTRIBUTION IS CLOSE TO 7 ~ 1 

Let us consider the electron energy distribution to be 



7e < 7i 
7e > 7* 



(20) 



We will take pi to be very close to 1 and P2 ~ 2. The exact value of pi can be found with i l 1 St . 
The synchrotron flux is given by 



1/3 



(21) 



(22) 



where A is a constant proportional to B, and 7„ D is the LF of electrons radiating at synchrotron frequency v , which is 
The frequency v(~{e) is the synchrotron frequency of electrons with LF 7 e 

= as-c (23) 

Substituting the last two expressions into l!21t . we can integrate this expression using l !20t . By absorbing constants in a new variable A', we 
find 



1/3 



Pi 



-pl + l/3 



P2 - Pi 

Pa - 3 



With this, we can find the spectral index for synchrotron radiation to be 



d In f Vo 
dlnu 



\ (Pi ^ I) 



P2-P1 Two ~\ 



\Pl-l/3 ' 



(24) 



(25) 



Since we are interested in the low energy spectral index, a, at u = 10 keV, 3^ = ( °^ J . For example, for v p = 1 MeV, then 

^ = 0.1, and for pi = (1.02, 1.1, 1.15, 1.3) and p 2 = 2 we find a = (-0.057, -0.089, -0.109, -0.173). Note that from ||25) a < 
forpi ^ 1. 

Since we want the spectrum at 10 keV to have a « 0, we will choose the LF of electrons to be y e = 74. We will determine how deep 
these electrons have to be in the KN regime (774) and their Compton-F parameter, Ykn, 4, so that we can obtain pi —¥ 1, We describe the 
calculation in the next section. 



5 PHYSICAL PARAMETERS CONSISTENT WITH a -> 

In this section, we determine the values of 7/4 and Ykn, 4, that are required in order to obtain the low energy spectral index close to zero at 10 
keV. The idea is simple: we scan all possible combinations of 774 and Ykjv,4 and determine the power-law index of the electron distribution, 
pi, which can be obtained with d 1 8b . Next, we use this power-law to find the observed spectrum at 10 keV, a, using i25\ . We present the 
results of our parameter search for 774 and Ykn,4 that yield a certain desired value for pi and a (Fig.[T] Left panel). 

For a given 774 and Ykn, 4 we can determine the radius, R, at which the emission is produced. Since Ykn, 4, is defined as Ykn, 4 = 
,kn* we can use an d {1} to determine R as a function of Ykn,4 and 774. We find 
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Figure 1. Left: Solid contour lines represent constant a (spectral slope at 10 keV), in the plane of the Compton-y parameter for electrons radi- 
ating at 10 keV, Y/fiV4> versus r)4, which indicates how deep these electrons are in the KN regime. From the upper right to the lower left corner, 
a = —0.1, —0.12, —0.15, —0.2, —0.3. The short-dashed lines are contour lines of the power-law index of electron energy distribution, pi, From the 
upper right to the lower left corner, p± = 1.1, 1.3, 1.5, 1.7. Right: Same contour solid lines of a in the Ykn,4 versus r?4 plane, but now the short-dashed lines 
represent contour lines of constant radius of emission. From left to right, log(_R) = 11, 12, 13, 14, 15, 16, 17 (in cm). The smallest radius is the photospheric 
radius, while the largest radius is approximately the deceleration radius (see text). Both figures are for the parameters of GRB080916C (details can be found 
in the next section). 



Ri 



3/2^2 r 1/2 

1 



Y, 



1/2 



(26) 



We can now calculate R as a function of r)4 and Ykna to find the range of radii where synchrotron radiation should be produced for each 
value of a. In the Ykn,a versus 774 plane, lines of constant radius are approximately given by Ykn,a oc ryf (Fig.[T| Right panel). 

It is important to mention that our analytical results use a simplified approximation for the KN cross section in l[8}, which is that 
<jkn ~ o't/t). However, results shown in all figures use the full KN cross section <jkn = ctt/(^), so that the radius is R oc r/ 2 (/(rj)) 1 ^ 2 , 
instead of R oc rf^ 2 ; analytical results are within a factor ~ 2 of numerical values. 

5.1 Constraint on the radius of emission 

The first constraint that we can place on the Ykn,a versus 774 parameter space is that the cooling time scale for electrons radiating at 10 keV 
cannot be larger than the dynamical time scale. After all, we need electrons to cool rapidly via IC scatterings in the KN regime to obtain 
a « 0. This constraint can be obtained in the following way. 

The electron cooling time in the jet comoving frame for electrons of LF 74 is given by t' cool — (l/t' 3yn + 'i/t' IC KN )~ 1 . We use l[6ll 
and l[8} to determine t' cool . We impose the constraint, as discussed above, that t' cool < t' dyn , and with it, we find an upper limit on the radius 
of emission, R coo i, and we can use it to constraint our Ykna versus 774 parameter space. An analytical estimate of this maximum radius is 
provided below. 

The cooling of electrons of LF 74 is dominated by the IC cooling in the KN regime. Therefore, we can approximate t' cool w t' IC KN 
and then set t' coo i < t' dyn , and use <fsj> and d 1 0b to find 



cool ,15 



220 



r 2 f ' 



(27) 



Another upper limit on the radius is given by the radius at which the external forward shock sets in, that is, the deceleration radius, Rdec- 
This radius is a function of the total blast wave energy, E, the circum-stellar density, n, which we assume is a constant, and the bulk LF (see, 
e.g., Sari, Piran & Narayan 1998) 



Rd 



130K 



l/3„-l/3 r -2/3 



The true upper limit will be given by the minimum of R coo i and Rdec- 

There is also a lower limit on the radius, given by the photospheric radius, 



Rph 



87rm p c 3 r 3 



(5.5 x 10 10 cm)L 5 3r3 3 . 



(28) 



(29) 



Therefore, the radius of emission should lie between R p h and the minimum of R coo i and Rdec- 
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5.2 Constraint on SSC component 

The last constraint that we place on the Ykna versus 774 parameter space is that the total luminosity in the synchrotron-self-Compton (SSC) 
component at ~ 1 GeV, Lie, should not exceed the synchrotron luminosity, L syn , for consistency with the Fermi data (Abdo et al. 2009), 
that is Lie /Lsyn < 1. The ratio of luminosities is given by Lie / L syn ~ Ykna, where Ykna is the Compton-Y parameter in the KN 
regime of electrons radiating at 7i (Nakar et al. 2009). The relationship between Ykna and Ykna, ls approximately given by 

V- V f'syn(< V KN (ji)) 

YKN,i = YkN,4-T, 7- ; 7 TV 5 (M) 

fsyn{< ^ K jv(74j) 

where f' 3yn (< ^A'jv(7e)) is the synchrotron flux below the KN frequency, v KN (-f e ) = m e c 2 /(h/je), both quantities in the jet comoving 
frame. The ratio of these comoving synchrotron fluxes is given by 74/7* = (10 keV/V p ) 1/2 , therefore 

r ra)i = y„, 4 M 1/2 . (3D 

We have found that the analytical calculation of Compton-F can overestimate its true value by up to a factor of ~ 10 (Barniol Duran 
& Kumar 2011, see, also, Nakar et al. 2009). If we restrict our parameter space to Lic/L syn < 1, to avoid conflict with Fermi high energy 
observations, then these last two considerations translate to a conservative constraint on Ykna given by Ykna < 10. Solutions that have 
Ykna larger than this limit violate Fermi observations and are ruled out. 

To summarize, we calculate the power-law index of the electron energy distribution function, pi, as a function of Ykn,4 and 774 using 
d 1 8 b - With pi and equation l !25t , we can determine the lower energy spectral index at 10 keV, a. The 2-D space (Ykna, V*) can be constrained 
by calculating the radius of emission that corresponds to each point in the Ykna—V* plane. The radius of emission should not be smaller than 
the photospheric radius nor larger than the minimum of the deceleration radius (eq. 128) . and the radius at which the dynamical and cooling 
time scales are equal, (eq.|27t. We can further constrain the parameter space by ensuring that the amount of energy in the SSC component is 
not excessive. 



6 APPLICATION TO GRB DATA 

We present the allowed Ykn,4,~V^ parameter space for two GRBs. We chose a very energetic Fermi GRB, GRB080916C, and another more 
"standard" GRB, which we call GRBt]. The parameters for GRB080916C are e p = 7.5, z = 4.3, L53 wl,T= 10 3 and, for most of 
the duration of the prompt emission of this GRB, a = —0.02 ± 0.02 (Abdo et al. 2009). To find the deceleration radius, Rdec, we choose 
E55 « 3 (Kumar & Barniol Duran 2010) and assume no = 1, although the dependence on E and n is weak. We choose the parameters for 
GRBt] to be more typical values of a GRB: e p — 3, z = 2, L53 ~ 10 -2 , T = 300, L55 rss 10~ 2 , no = 1 and, for most of the duration of 
the prompt emission of this GRB, a = (Preece et al. 2000). Following the prescription at the end of last section, we present the results of 
the allowed parameter space in Fig.[2](Left panel). Note that we present the parameter space now as a function of Ykna, instead of Ykna 
(which is a better indicator of Lie / L svn ), however, they are related by a constant factor, eq. d3H . 

As can been seen, a —¥ as we move into a region where 774 and Ykna (and consequently, Ykna) both become larger. What is the 
maximum possible lower energy spectral index at 10 keV (a rnax ) that can obtained? It is given by the value where the maximum of Ykna 
intersects with the maximum allowed radius (see Fig. [2] Left panel). This intersection gives us the maximum allowed value for a (a max ) 
which is found to be —0.1 (Table QJ. Therefore, any of the GRB gamma-ray radiation for which a > —0.1 cannot be produced by the 
synchrotron process in shock heated plasma, where electrons are only accelerated when they cross the shock front and are scattered back to 
the other side (this however, is not the case if electrons are continuously accelerated while they are traveling downstream or upstream; we do 
not consider this scenario in this paper). 

In Table [T] we also present the maximum value of a obtained if: (i) We decrease the value of the peak energy as observed during the 
prompt phase (Abdo et al. 2009), and/or (ii) We decrease the value of the LF of the source, as suggested by several groups (see, e.g., Zou et 
al. 2011, Hascoet et al. 2011). 

The luminosity carried by the magnetic field, Ls, as measured by a lab frame observer is 

L B = ?^—4ivR 2 c = (1.1 x 10 52 erg)R 2 5 (e p /n 4 ) 4 r3 4 (1 + z) 2 , (32) 
Sir 

where we made use of (|5j- Therefore, the fraction of energy carried by the magnetic field is 

l b A1 , , >4 r i (Ri 5 {l + z)\ 2 
eB = — = 0.1(ep/7H) L 5S I ^2 I . (33) 

We present the value of e_g for the a rnax case (TableQ3; es is found to be very small, in the range ~ 10~ 6 - 10~ 4 . We also provide in Table 
1 the value of the LF of electrons that radiate at the peak of the spectrum (7^) in order for the low energy spectral slope to be a max ; we find 
7* > 10 6 . 

If we take the observer frame variability time scale of the gamma-ray light curve to be 
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Figure 2. Left: The solid lines are contours of constant a (from the upper right to the lower left corner a = —0.1, —0.12, —0.15, —0.2, —0.3), while the 
short-dashed lines correspond to constant R value in this 2-D plane; R = 10 11 cm for the left most line, and the subsequent ones correspond to 10 12 , 10 13 , 
10 14 , 10 15 , 10 16 and 10 17 cm (the largest radius is approximately the deceleration radius). Note that we plot Yk N,i instead of Yjc jv,4 (which differ only by 
a constant factor). The maximum allowed value of Ykn i ~ 10* so as to avoid a SSC component excess in the Fermi high energy observations, is plotted as a 
long-dashed line (the factor of 10 allows for an overestimation in the analytical result). The allowed region is the region below the long-dashed line and to the 
left of the maximum radius, which has been shaded. The maximum value of a max = —0.11 is obtained at the intersection of the maximum allowed radius 
and the maximum Ykn,i- It is clear that a > — 0.1 is impossible to reach for the data of GRB080916C. Right: The solid lines show a in the e f /eg versus 
r}i plane. e e and eg are the fractions of energy in electrons and magnetic field, respectively, while rn indicates how deep in the KN regime are the electrons 
radiating at the peak frequency. From top to bottom, a = -0.20, -0.25, -0.30, -0.35, -0.40, -0.45. The data in these panels are for GRB080916C, 
however, the plot for an average long-duration GRB is almost identical. 
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Table 1. We present the maximum possible spectral slope, a max , for two GRBs: GRB080916C and GRBi]. a max is the value of a obtained, where the 
maximum Y^N,i intersects with the maximum radius (see Fig. [2] Left panel). We present results for different values of the GRB-jet LF and the observed 
spectral peak (e p ; see eq. 1 for definition). In addition, we also present the values for rjn, 7^ and eg, where the spectral slope is maximum. The value of a does 
not reach the observed value a = —0.02 ± 0.02 (Abdo et al. 2009) for GRB080916C, nor a = for GRBt] for any combination of parameters. At a max 
we find extreme values for the LF of electrons radiating at the peak and for the energy fraction in the magnetic field. The values of ct max where we assume a 
variability time scale of St = 0.1 s are in parenthesis (see next Section). 



as expected in the internal shock model, then in order for the low energy spectral slope to be a max , we find St ^ 10 s, since the emission is 
produced at (or very close to) the deceleration radius. This is very long compared with the observed time scale of 0.1 s or less. 

We can calculate the low energy spectral index as a function of e e /es, where e e is the fraction of energy in electrons. This will allow us 
to compare our results with previous work that use this ratio. To do this, we use (Ando et al. 2008) 



Y KN ,i = \l^^-, (35) 
rji 

where r/i parameterizes how deep electrons of 74 are in the KN regime, analogous to $3$, and it is related to 774 as ?7;/?74 = 71/74 = (v p /10 
keV) 1/2 . Using (3D and {35J, we find e e /e s as a function of 774 and Ykn,a and calculate pi using equation d!8t . and a from equation d25t . 
The result is shown in Fig. |2] (Right panel). 
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7 DISCUSSION 

Most GRBs have low energy spectral index a between and —0.1; f v oc v a . For a > —0.1, it is required that the power-law index for 
electron energy distribution function, pi = | ^j°" e |, should be less than 1.15 at 7 e corresponding to 10 keV synchrotron photons (in observer 
frame). For this, two conditions must be satisfied. 1. The Compton-y parameter for electrons radiating at 10 keV should be Ykn,a > 20 
(including KN effect), and 2. The LF of electrons radiating at 10 keV should be > 2000re~\ where e v = ^^JT' - , that is, T]a > 2000 (see 
Fig.[TJ. These conditions apply both for the data of GRB080916C (a highly energetic explosion) and also for an average long duration burst. 

Electron index changes sharply, from being close to 1 (when IC in the KN regime dominates) at 10 keV to > 2 at the peak of the observed 
spectrum at ~ 1 MeV, that is, the electrons distribution index increases from ~1 to >2, when the electron LF increases by a factor ~ 10. In 
this case, the spectral index for synchrotron radiation at 10 keV is not given by —(pi — l)/2 (as mentioned in the Introduction); in fact, it is 
significantly smaller. For instance, when pi = 1.02 (for electrons radiating at 10 keV), = —0.057 and not —0.01 as naively expected; 

see ( 1251 ). 

Another consequence of the requirement that r)4 > 2000 is the that energy fraction in magnetic field, es, is rather small es < 10 _9 i?f 5 
(GRB080916C) and eg < 10 _7 i?i 5 (for an average burst); see ( 1331 , In addition, the LF of electrons radiating at the peak of the spectrum 
(~ 1 MeV) is ~ 10 times larger than that of electrons radiating at 10 keV, therefore, ji > 2 x 10 4 re~ 1 , which is ji > 10 6 almost independent 
of GRB energy (Table 1). Both values of es and ji are extreme and their implications will be discussed in the Conclusions. 

However, we can also ask: what is a for a reasonable set of parameters? In this case, "reasonable" means two things: (1) The radius 
of emission of the prompt emission should be between (i) the photospheric radius and (ii) the radius where electrons producing 10 keV 
synchrotron photons cool on a time scale shorter than the dynamical time via IC scatterings in the KN regime, or the deceleration radius, 
whichever is smaller, and (2) The energy in the IC component should not be very large, so that the GRB spectrum does not show an IC 
bump at ~ 1 GeV, as the Fermi satellite sees no sign for such an excess. We have calculated the maximum value of a that can be obtained 
for a highly energetic Fermi burst, GRB080916C, and also for an average GRB and the results are presented in Table Q] We find that the 
maximum value of a is a = —0.1. For this value, the same consequences as discussed above apply, and can be found in Table Q] Namely, 
that r/4 ~ 10 3 , and this implies a very large value of -ji and an extremely small value of es. Moreover, the maximum value of a occurs at 
a very large radius, close to the deceleration radius, which, will have problems producing variable light curves with St smaller than a few 
seconds. 

Conversely, we can also fix the observed variability time scale of the gamma- ray light curve to be St = 0.1 s and determine the radius 
of emission with d34t . as expected in the internal shock model. We can determine the maximum value of a at this radius and its value is 
presented in parenthesis in Table [T] Notice that this value of a is even further away from the observed value. However, at this radius, the 
values of e B ~ 10~ 4 - 10" 3 and 7, - 10 3 - 10 5 are less extreme and in marginal agreement with the internal shock model. Nevertheless, 
this further accentuates the fact that synchrotron emission, in which electrons cool mainly via IC in the KN regime, in the context of the 
internal shock model, cannot explain the observed a ~ spectrum of most GRBs. 

Expressing our analytical results as a function of e e /es allows us to compare them with previous numerical work. Nakar et al. (2009) 
(see, also, Fan 2010) have found numerically that with e e /es = 100 (10 4 ), a cannot exceed ~ —0.3 (—0.2). In this work, we analytically 
confirm their results (see right panel of Fig. |2). Daigne et al. (2011) have also found numerically that a = is possible in the case when 
electrons are injected instantaneously. However, in Nakar et al. (2009) and the present work, we have assumed that electrons are injected 
regularly over a time scale comparable to the dynamical time scale. Nevertheless, Daigne et al. (2011) present numerical results for our 
scenario in their fig. 2 (bottom left panel) for r\i — 100, which agree with our analytical calculation. 

Our work and the work by Daigne et al. (2011), as mentioned before, differ mainly on the chosen time scale at which electrons are 
injected to the shock, ti n j ec . We take tinjec ~ tdyn, where tdyn is the dynamical time scale, whereas Daigne et al. take tinjec >C tdyn- We 
find a stricter limit on the allowed value of a: a < — 0.2, whereas Daigne et al. find that solutions with a — can be reached (see their fig. 
2). These results are not in contradiction, since both studies probe two different phases found in the internal shock model. The first phase 
(tinjec ~ tdyn) has electrons continuously being injected to the shock as it crosses the shell. The second phase (tinjec <C tdyn) corresponds 
to the case where the shock has already traversed the shell and electrons cool as the shell adiabatically expands. Most of the available energy 
is dissipated in the first phase; for this reason we have chosen this particular scenario in this paper, which is the common practice in studies 
of the internal shock model (see, e.g., Piran 1999). 

However, when calculating synthetic GRB light curves, Daigne et al. (201 1) do follow the dynamics of the shock crossing numerically, 
and consider a large number of discretized shells on the dynamical time scale. In each collision the electron injection occurs instantaneously, 
tinjec <C tdyn, but the electron injection process over the full simulation is comparable to the dynamical timescale. In this case, when 
including IC cooling in the KN regime (see their fig. 9), they find a < — 0.1. The difference in our results appears because the IC scatterings 
in their work do not occur between the same photon and electron distributions we have considered: In Daigne et al. (2011) the scatterings 
between photons emitted in a shocked region and electrons or photons present in a subsequent shocked region were not considered (see 
Bosnjak et al. 2009). This affects the cooling of the electrons via IC, allowing them to reach spectra closer to a = 0. 
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8 CONCLUSIONS 

In this work we have investigated the possibility that the observed low energy GRB prompt spectrum, which is f u ~ i>° (below the peak) for 
a good fraction of all long duration GRBs (Preece et al. 2000; Kaneko et al. 2006; Pelangeon et al. 2008; Krimm et al. 2009; Ghirlanda et al. 
2010), is due to synchrotron radiation from electrons that cool mainly via the IC mechanism in the KN regime (Derishev et al. 2003, Bosnjak 
et al. 2009, Nakar et al. 2009, Wang et al. 2009, Fan 2010, Daigne et al. 201 1). 

We present an analytical method to determine the power-law index of the electron energy distribution function, pi, that cools via IC 
cooling in the KN regime as a function of two parameters: 77, which is a measure of how deep electrons of interest are in the KN regime, 
and Ykn, which is the Compton-F parameter (including KN corrections) for these electrons. We have calculated the observed low energy 
spectral index for synchrotron radiation, a, as a function of these two parameters as well as the power-law index of the electron energy 
distribution above 7s (7^ corresponds to the LF of electrons radiating at the peak). We find that a is not simply given by — (pi — l)/2 as 
naively expected, but it is smaller, which makes it very difficult to explain the observed value of a ~ for a good fraction of GRBs. 

We find that a > —0.1 cannot be obtained for parameters relevant for GRBs, if the radiation mechanism is the synchrotron process and 
electrons are accelerated in a relativistic shock, where electrons are only accelerated when they cross the shock front and are scattered back 
to the other side. Therefore, the 7-ray radiation from a significant fraction of long duration GRBs that have low energy spectral index larger 
than —0.1 cannot be accounted for by this mechanism. 

Even q ~ —0.1 faces severe difficulties. The large radius for generation of 7-rays is in conflict with the short variability time ( < 0.1 s) 
of prompt GRB light curve. Moreover, the energy in the magnetic field must be extremely small, es ~ 10~ 6 — 10 -4 , and ji ^ 10 6 for the 
mechanism to be able to able to harden the spectral slope from a = —0.5 to ~ —0.1 (Table 1). 

It is unlikely that the energy fraction in the magnetic field will be so small (es < 10~ 4 ) in internal shocks. If the central engine of GRBs 
is powered by accretion onto a black hole, we expect es ~ 1% as magnetic fields of such a strength are likely produced in the accretion 
disk by the Balbus-Hawley mechanism (Hawley, Gammie & Balbus 1996); for a magnetar based central engine this small es is even more 
surprising. 

For the typical thermal LF of electrons to be large, 7, > 10 6 , in internal shocks where shells collide with a relative LF of a few to 10, 
it is required that approximately only 1 in ~ 10 3 electrons are accelerated when they cross the shock-front but they receive ~ 10% of the 
total energy. This is in contradiction with the numerical PIC simulations of Sironi & Spitkovsky (201 1). Moreover, the ~ 99.9% of electrons 
which are not accelerated have a thermal LF of a few thousand due to their interaction with protons (Sironi & Spitkovsky 201 1), and these 
electrons produce a significant IC bump in the spectrum at ~ 100 MeV which is not seen for any bursts. The SSC flux of these electrons at 
100 MeV will be very large: about a factor of 10 larger than the observed flux. 

All these difficulties suggest that the synchrotron radiation mechanism in internal shocks does not provide a self-consistent solution 
when the low-energy spectral index for GRBs is larger than about —0.2. 
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